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Abstract 



<J->. 

For the large-angles hard photon emission by initial leptons in process of high energy 
annihilation of e + e~ —> to hadrons the Dirac tensor is obtained, taking into account 
the lowest order radiative corrections. The case of large-angles emission of two 
hard photons by initial leptons is considered. This result is being completed by 
the kinematics case of collinear hard photons emission as well as soft virtual and 
real photons and can be used for construction of Monte-Carlo generators. 
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1 Introduction 



The problem of precise knowledge of the cross section of annihilation e + e~ to 
hadrons caused by the long staying problem of theoretical estimation of muon 
anomalous magnetic moment g — 2 [1] : 



DO 



„ 3\ 



4m^. 

a e-e^adr X 2 (l - x) S 

R(s) = — z — — , K '(s) = / ax— z ; -, p = — zr- 

Extraction of cross-section e + e~ — > 7* — > hadrons from experimental data is 
one of the main problem of modern experimental physics. The Monte-Carlo 
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programs creation which takes into account the emission of real photons by 
the initial leptons is the motivation of this paper. 

Dirac tensor (cross-symmetry partner of Compton tensor) i.e. the bilinear 
combination of the currents of hard photon emission averaged on leptons spin 
states and summed on photon polarization states takes the contribution on 
Born level and the ones arising from 1-loop correction. Infrared divergences 
are parametrized by the introduction "photon mass" A. In the final answer 
it is removed in a usual way by adding the contribution from additional soft 
photon emission. 

We don't consider photon emission by the final charged particles as well as 
the effects of charge-add interference of emission of virtual or real photon 
emission from leptons and hadrons. So the Dirac tensor obtained in such way 
is universal. 

The paper is organized as follow. In the part 2 the relation of Dirac tensor with 
cross section of the radiative annihilation of lepton pair to hadros is clarified. 
We put the Born level expression for Dirac tensor and rerive the general form 
of radiative correction to it by using the symmetry relation. 

In part 3 we obtain the contribution arising from mass operator of positron 
and vertex function for the case when positron and photon are on mass shell. 
In the part 4 we consider the contribution from vertex function for the case 
of electron on mass shell and the box-type Feynman amplitude with electron, 
positron and one of photons on mass shell. 

In section 5 we analyze the total result for Dirac tensor, adding the emission of 
additional soft photon contributions, which provide the infrared divergences 
free final result. Limiting case of almost collinear hard photon emission is 
considered and some numerical estimation are given. 

We put the form of hadronic tensor for several final states: 7* — > 7r + 7r~, 
p + p~. In Appendixes A and B the details of calculation are presented. In 
Appendix C the contribution to Dirac tensor for the case of two hard photon 
emission is given. 



2 General analysis 

The Born level matrix element of hard photon emission by initial leptons in 
process of annihilation e + , e~ to hadrons through the single virtual photon 
intermediate state 

e + ( P+ ) + e"(p_) ->• 7*(g) + 1 (p 1 ) 1 (p 1 ) + %) (2) 
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has a form (see Fig. 1) 




(a) 

Fig. 1. Diagrams contributing in Born level. 



M = (±^- {p+)0( B) u{p )Hp{ql (3) 
n (B) P- - Pi « . ~-P+ +Pl 

p -x- -x+ 

where e(p\) is polarization vector of the real photon. The H p (q) is the current 
describing the conversation of virtual photon whith momentum q to hadronic 
state. We will restrict ourselves by kinematics conditions of large-angles scat- 
tering. 

s = 2p + p_, x± = 2 PiP±, Pi = °> P± = m2 , ( 4 ) 
s ~ X+ ~ X- = <1 2 , Q 2 > 0, s ~ q 2 ~ x+ ~ X- > ™ 2 - 

In expressions below we put m = everywhere except the denominators of 
loops integrals. 

Cross section can be expressed in terms of the summed on spin states of the 
module of matrix element square: 

E W? = (4™) 3 ^f^, 

spin \" / 

B ppi = ±Tr p + O p p_O pi , H ppi = ]T H p (q)H* pi (q). 

spin 

The differential cross-section can be written as: 



dr / = (2^(p ++ p_-p 1 -E^)n^ p - (5) 



For differential hard photon cross section we obtain 



Wider 6 e ^ x 2a 3 



d 3 pidrj s(g 2 
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where 



B PPl = B g g PPl +B^ +p p- n +B-^- p p- Pl +B + 4p-p + ) ppl , 

(p+P-)opi = P+pP-pi + P+piP-p- (7) 

The quantities with the "tilde" are defined as 

g PP1 = g PP1 - -^QpQpi, p±p = p±p - ip-^r- ( 8 ) 

In the Born approximation (see Fig. 1) we have 



B? = — (V + )d + X 2 J, B* + = B B _ = ^, B* = 0. (9) 
X+X- X+X- 



For q 2 = we reproduce the Dirac cross-section of e e + — > 77: 



dr 2a 2 xl + X 2 - 



dCi s x+X- 



(10) 



Below we concentrate on calculation of one-loop radiative correction to Dirac 
tensor. 

Let us now show that in considering of the corrections one can restrict only 
by half of full set of Feynman diagrams for considering process (2). So we put 
O p = + 0+ separating the contribution of emission from electron leg 0~ 
and positron one 0+ (see Fig. 1 for the Born case and Fig. 2 for the 1-loop 
corrections). 

One can show that using the cyclic property of the trace as well as the mirror 
property 

Tr aia 2 • • • a 2n = Tr a 2n • • • a 2 ai , 
that the total contribution to the Dirac leptonic tensor can be written as: 



= (1 + A m )(l + V)Trp + 0+p-6 B pi . (11) 

Here the exchange operations acting as 

A F = F 
^ppi 1 PP i 1 p\pi 

VF(p + ,p^,p 1 ) = F(-p_, -p + , -p^, 
VF(s, q 2 , X +, X-) = F(s, q 2 , X -, X+) = F. (12) 

Here and below we imply only the real part of leptonic tensor. 
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3 One-loop corrections. Real photon vertex and self energy con- 
tribution 



The virtual correction of lowest order is described by 8 Feynman diagrams 
shown on the Fig. 2. 








(e) (f) (g) (h) 

Fig. 2. Diagrams contributing in 1-loop level. 
Let us distinguish contribution of FD Fig. (2(e-h)) to 3 classes 

Trp + o+ P _d£ = + t;;- + rf pi 

with T box and T vcrx correspond to Fig. (2 (h,g)) and T s to Fig. (2 e,f). 



(13) 



Consider first the contribution to matrix element arising from Feynman dia- 
gram Fig. 2 e,f. 

Matrix element of FD Fig. (2 (e)) contain the mass operator of electron 
In kinematics conditions of our problem (x+ ^ m2 ) we obtain [3]: 



a ,3 1 



2tt v 2 2 



X+ 



(14) 



/ 1 *± 7 1 m 



where A is the so-called " photon mass" . 

Matrix elements of FD Fig. (2 (f)) contain the vertex function with real pho- 
ton [3]. 



M 



' = s-H/ mm — u[p - y 



We use here the notations (see eq. (A.l)) 
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Using the relevant loop integrals, obtained in [2] (see Appendix) we have a 
matrix elements of FD Fig. (2(f)) which contain the vertex function with real 
photon [3] 

M f = ^v(P + )[- — (h ~ hpie + e(l x - \l + - h] ~ P++Pl l P u(P-)- (15) 
^ X+ A z z ~x+ 



As a result we obtain the infra-red free and gauge- invariant expression: 

a 11 

M e + M f = -^ + v{p + )p 1 e lp u{p^), $ + = — - Z+ - - . (16) 
71" 2x+ 2 



Inserting this expression to the relevant part of 0+ we obtain for 

Tp pi = -$+Trp +: pi 7A7o p_ [— 7a (p_ - pi)^ + — lpi (-p + + p x ) 7A ] , (17) 

X— X+ 

where we used relation p ip = {p + +p~) p keeping mind the gauge invariance of 
hadronic tensor q^H^ = 0. Expression (17) leads to the form for contributions 
of diagramms Fig. (2(e,f)): 



2 

Tf pi = -^ + [2p_ p p_ pi (^- - 1) + 2p_ pP+pi (^- - 1) - (s - X -)g PPl ]- (18) 

X — X — 



Applying the operation 1 + A ppi and 1 + Pwe obtain the full result: 

(l + A m )x(l + P)T p s pi 

= Bgg ppi + + B^p +P p +P1 + Bf_ ) PP1 , (19) 

with 

Bf = ^[sc - Xl ~ Xl% + T= B- = - X _]l s + T_ s , 

5? = ~—k 2 - X + ]*. + = [g 2 + s]l s + T + s _, 

X-X+ X-X+ 

c = X + + X- (20) 

Here we use the notation 

Tf = -2^^[1 + 2l sp ] - 2^*±[1 + 2U, 

X+ X- 

T S = -^—[q 2 - + 2/ sp ] , 7? = — ^" [? 2 - X+][l + 2 U; 

X+X- X+X- 
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tL = -A_ [ s - x _][i + 2 y + -4-[ s - x+][i + 2i s 



X+X- 



X+X- 



sp 



If, — L 



h - L 



(21) 



4 Vertex and box type diagram contributions 



Contribution of Fig. (2(g,h)) can be written as 

1 ppi ' rl ppi — , . ..2 ' 



X- X+ X~X+ X+ 



with 



Si-- 
S 2 ~- 



1 d 4 /c Tr B p p-j v {p- - £1)7, 



pi 



4» 2 (0)(2)(2)(g) 
ld 4 /cTr J B p p_7 Pl (-p + + pi)^ r , 
4^ ^ (0)(2)(2)(g) ' : 
1 d 4 /c Tr V p p_-f v (p_ - pi)7 pi 



C, 



rid 
~ J 1~i 



4i** (0)(2)(g) 
1 d 4 k Tr V p P-^f pi + £1)7,, 



(0)(2)(?) 

B P = p+ix(-p+ - k)l v (-p++Pi - k)l P (p- - k)l\, 

V P = P+lv(-P+ + Pihx(-P+ + Pi - k)l P {P- ~ (23) 

Using the loop integrals listed in Appendix A we obtain: 



T iT + T£ = D g g PPl + D_p_ p p_ pi + /J+p+pP+p, 

+ D+~P+ P P-pi + D~+P~ P P+pi , (24) 

Applying the interchange operator 1+A ppi and 1+V, -D(x+, X-) = ^-Kx-j X+)> 
and rearranging the gauge-invariance we put it in the form: 

(1 + Ap Pl )(l + V)(T^ + T p v p f ) = B™~ 9ppi + BVf-p-n 
+Bl B p +p p +Pl + Bl B (p + p_) PPl , 

where 

B v g B = 2{D g + D g ), B VB = 2(D_ + D + ), B v + B = 2(D + + L>_) 

= + + /3 + _ + £_+. (25) 
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Here we could see by construction that 



B v g B = B v g B , B VB = Bl B , Bl B = B VB , (26) 



and in explicit form B\ are 



= 4 ££ -8 f!;s + 2*1+2x1-4.0 + 4,' , + _ + ™ 



= 8 -^^l s + -^[ll + 2(/ s - 1)/ A - Q + TY B , 



X-X+ x~x+ 

B v + B = 8 -^^l s + -^[l] + 2{l s - l)l x - l s ] + T™, 

X-X+ X-X+ 



B\1 = 4(s + g \ + T™, 



X-X+ 



where the expressions T^ B contains non-leading terms. 

These quantities contain the ultra-violet cut off logarithm L = ln^ which is 
eliminated by standard regularization procedure [3] L — > 21 \ — 9/2. 

Collecting the leading terms whcih contains the large logarithm l s and infrared 
one l\ we obtain: 

{B^ B + B^)i eading = 2B B (l 2 + 2(l s — l)L\ — 3/ s ), 
(B+ B + B^_)i eading = (BY. B + B^)i eading = 2B B (ls 2 + 2(l s — l)L x — 3/ s ), 

(B+ B + B+_)i eading = 0. (27) 



5 Discussion, explicit form of tensor structures. 



The infrared divergences constrained in contribution of virtual photon emis- 
sion canceled when takes into account the emission of additional soft photon 
(center-of mass of e + e~ initial is implied) 



d(T Lft = $soft<icr B , 

= - is? / + & • w < A£ ^ A (28) 

where w = Vk 2 + A 2 . Using the standard integrals we obtain 

<W = - 1)('a + 21n^) + i/ s 2 - y]. (29) 
Summing all contributions we find Dirac tensor: 
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B ppi = ( B g 9ppi+B + ^ + pP~ Pl +B_J- p p- Pl ) 

a., w 3 , AE . a. n 2 3.. 
x (1 + -{I, - 1 - + 21n— + -(-- + - 

7T Z £/ 7T O 2 

Oi 

- [Tg9p Pl + T^ p p^ Pl + T + p +p p +Pl + T+_ . (30) 

Quantities = Tf + T^ B are free from infrared singularities and do not 
contain large logarithms. Quantities Tf are given in (21). Quantities Tf B are 
given in Appendix B. 

Expressions for Tj contains in addition nonphysical singularities x± 2 j X± 3 - 
Nevertheless one can be convinced in cancelation of terms proportional x± 3 
with structure G (see eq. (B.2)) and terms x± i n convolution of Tj. To be 
defined let us consider the case of small values of x~, m 2 <^ X- ^ s ~ Q 2 ~ X+- 
It corresponds to the kinematics p\ = yp-. In this case non- leading terms 
containing poles could be put in the form: 

T g g ppi +(T- + fT + -2yT + _)p_ p p_ pi , y = l-y, y = *± (31) 

s 

and this combination contains only the lowest order pole X- 1 - The Dirac tensor 
in the limit m 2 <C X- ^ s has a form: 

B% = B^(l + ^(l s - l)( 3 - + 21n^) + ^ - £)) 



where 



Kpl im=l - ± ^(9 P p, + -P-pP-p 1 ), x = l-x, x X 



)B, Urn _ 

xy "" ' S " ' " ' " S 



pl ° m = v y ( 16 " 18y + 1%2 + 4(1 + v2) [ lnVln l - Li2 (rh;) _ t 



-8yylny + (-12 + 20y - Uy 2 )\ny , (33) 



Values of xyT l ^ m as a function of y and x = 0.1 are presented at Fig. 3. 

Obtained formula have a power accuracy as well as we systematically omit 
the terms of order m 2 / s compared to ones of order of unity. 

Similar properties have a cross-channel-Compton tensor with one real and 
another virtual (space-like) photon [4], where terms of order m 2 / s was taken 
into account. 
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Fig. 3. xyTg as a function of y at x = 0.1 
6 Samples of hadronic tensors. 

Hadronic tensor is the summed on spin states of bilinear combination of matrix 
elements M p M* p ^ , where the current M p describes the conversion of heavy time- 
like photon to some set of hadrons. For the case of creation of a pair of charged 
pseudoscalar mesons (it + tt~ , K + K~ , ...) we have 

H p P k pt = (p s + - p s jm - P-)pi ,q=p s ++p s - 

For conversion to pair of charged spin 1/2 fermions 7 — > p + (p+) + p~{p-) we 
have 

H%T = 4K> m pi + P^P- P - j9 PPl }. (34) 
For creation of a pair of charged vectors mesons p + p~, K* + K*~ one obtaines 

(35) 
(36) 

The gauge invariance requirement H ppi q p = H ppi q p = is fulfilled. 



H*g- Wig, 



q 2 



+ (q+ - q-) P (q+ - g-) w (3 
,2 



T 

m p 
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A One-loop Feynman integrals 

In this section we perform the result of calculation of 4-fold integrals, asso- 
ciated with one-loop Feynman diagram. Here and below we imply only real 
part of integrals. The denominators of integrals defined as 



(0) 
(2) 
(2) 
(?) 



k 2 - A 2 , 

(p_ - k) 2 - m 2 + iO = k 2 - 2p_k + iO, 
(-p+ - k) 2 - m 2 + iO, 
(pi — P+ — k) 2 — m 2 + iO. 



(A.l) 



The four denominator scalar integral 




dk 



17T 



2 



(A.2) 
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has the form 



-^0225 — 



SX+ 



L 2 - 2l + l s - l„k + 2Li 2 [ 1 - j - 



(A.3) 



where the logarithms was denoted in (2) and 



l q = ln-^r, l s = \n^—. 



(A.4) 



For the tree and two denominator scalar integrals 



r 



(A.5) 



with r = (ij), (ijk), (ijkl) where k, I = (0), (2), (2), (q), we have following 
expressions 



1 



l 02q 



L 22q 



L 02q 



2Xh 



2tt 



2i 



1 



2(s - g 2 
1 



3 

H - H 



I022 — 77" 
2s 



/ s 2 + 2l a l t 



4tt 2 



Z 9 (^-Z + ) + ^(^-M 2 + 2Li 2 fl + ^ 



3tt 2 



X+ + Q 2 

Two denominator scalar integrals are 

/,„ = L - In + 1 



I 2=L + 1, 
J o2 = L + 1, 



2g - ^ - „ ? 

I 2 - 2 = L-L S + 1, 



(A.6) 



J 0g 

4, = i " 1- 



The vector integrals can be defined as 

d 4 kk" 



I 



a+ql + a r qt + ajpj*. 



(A.7) 



For the vector integrals with two denominators we have (imaginary part is 
neglected) : 



a 2q - a 2q ~ ~ a tq ~ 2^ ~ ^ + 2 J 



"22 



~ a 22 — o( L - h + 7; i- 



~ - l T 1 

a 02 — 4' 



a 9 = 


-4q = 




a 2, = 


-\4q - 


i 

" 2 






1 


u 02 


4- 


4 



L-L + 



27' 



(A.8) 
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and the coefficients for the vector integrals with three denominators are 



a 02q — 
a 02q = 

a + - = 

U 02q 

a 22q = 
a 22q = 



- X+^02g + 1+ + lq , 

a\ a a 



a 02q ~ a 02q ~ _ [ W lq ) ) 



l 02< ? 



~~ I()2q l+i 



c lg ) > 
s r 1 

_ -<22g H 

C C 



a = x+ + 9 , 



°022 °022 



'1st 



"22q 



~^22q + - ( h — lq ) , 



2s 



l q + 2)-—(l a -l q ), c = s-q 2 = X+ + X- (A.9) 



Finally, the coefficient of the vector integral with 4 denominators has the form 

a 1 = ^ (x+A + X-B -sCj, a+ = T { X+A " X - B + S °) 

a~ = ^(^-x+A + X-B + sC^j, d = 2s X+ X-, 

B = I()2q — hlqi 



A — l22q — I()2q, 
C = lQ2q — 1 022 — X+Io22q- 



(A.10) 



The second rank tensor integrals can be parameterized in the form 
huh,. 



T\iV 

r 



a 9 r g + a^piPi + a+ + g+g+ + a r 1-1- + a l + (Pi1+ + Q+Pi) 



. (A.11) 



\1V 



The coefficients for tensor integral with four denominators are (we suppressed 
the index 022g) 



- (a 6 + A 7 - A 10 ) , 

x+y ) 

a 1 - = -j- (A 2 + A 7 - A w ^j , 
A 5 - X+a 1 ^, 



X- 
1 



a = - \A 2 + A 6 - A 10 

s 

1 



a L = — [Ax - sa 



X- 
1 



= - s (A 3 - X -a 1+ 



A 10 - A 2 - x+a 



i+ 



with 



A, 


~ a 22q a 02qi 


A 6 


= a 02q 


U 22g> 


A 2 


= a 22q-> 


A 7 


= a 02q 


- X+a 1 , 


As 




A 8 


= a 02q 


— °022 ~~ X+ a 



(A.12) 
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A 4 
A, 



a 02q a 22qi 
a 02g ~ a 22q-> 



A, 
A 
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l 02q ' 



a 022 X+ a 



(A.13) 



For the tensor integrals with three denominators Iq% q we have coefficient^!] 



a 02, 



4a 



l 02q 



-a 



029 ~ 2a 



(z + -g-i 



a 02g — a 02q ~ a 02q ~ 2a^ Q ^ 



L 02q 



2 T Jxl, (q 2 ) 2 + 4g 2 x + - 3x1 

X+ 1 02q H 1+ 

a 



2a 



g 2 + 3xh 



(A.14) 



The coefficients entering into the tensor integral I/tg are 



a 022 ~ TV^ — + 



-, a Q22 a 022 2s ^ s a ° 22 



1 

2? 



(A.15) 



and the coefficients for the tensor integral I^ q are 



"02q 



%2q 



Ah 



u- 



l 02q 



2*- 



■H++2), 



^029 - J 02< ? 



/o2 9 + ^(3i + -l). 



2X 



(A.16) 



In the case of the tensor integral J 22(? they have the form 



a 



22q 

+_+ 
22q 

h 22q 

b 22q 

,11 
L 22q 



1 r 3 S i Q i 

-L H i a + —L 

2 4 2c 2c 9 



"22g 



1 S S 

2 + Tc 3 ~ 2c q . 



— ~ S ^22q + —I 



5s. 2q 2 + 3s 



4s - q 2 + s 2 /^ - — Z s + 



2c 



2c 



2c (lq h 



22 i 2c 



3s 2 3s 2 - (q 2 ) 2 + Asq 2 



2c 



(A.17) 



(A.18) 



(A.19) 



Formula (2.15) form [2] contains a misprint 
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B Explicit form of coefficients of non-leading tensor structures 



T, 



J B — (1 + V) a + CLlhq + a 2^qp + a dsp + ^sqhp + 

+a 9 Li 2 1 



(B.l) 



where 



a = 



ai 



a 2 



04 



a 6 = 



a 8 



a 9 



7T 

T 

4s 



■10— + 22— 



X- 

+ io^± - 20 



c X- 
6(s-X+) 



X- 
s 

— + : 

X- X+X- 



X+X- 
8s 2 



<? 2 + x+ 






4x- 4x+ 


4s 
+ 


4s 
f — 


X+ X- 


X+ 


X- 


4y_ 8s 
— ^_ + 


4s 


8s 2 


X+ X+ 


X- 


X+X- 



2X 



+ 



X- 



4s 2 



6s 

— + 

X- X+X- 
4s 8s 8s 2 



aio = 4 + 



G = 



X- 
Li 2 



X^ 



X- 
+ Li 2 



X+X- 



+ 



,x+ 

X- 



K+ + 

X- 



1 ~l~ 9 ] "I" ^sq^sp 
Q J 



X+X- 



-z 2 + - 

2 S9+ 6 • 



(B.2) 



Note that in the limit X- ~> quantity G turns to zero. 



yVB _ jyrpVB 

T VB = bo + 6i/ sg + + Msm + Map + b^l sq l sm + b§l sq l S p + 67/^ 

+ 6 8 Li 2 ( X " 7) + 6 9 Li 2 (l + ^p) + & ioLi 2 + |r) 

- 8 (S ~ X + )3 G + s4 fl - — 1 PG, (B.3) 

X+X- X+ V X-/ 
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where 
b = 
+ 
+ 



c V x+ x+J 
2 r 44 44 4 X 

— + — + -v 

X+ X- X- 



A'4 



Q 



-8 + 



fXj 
X- 



+ 



32 

X+ 

-44- 

4s 

X+, 



48 

X- 



+ 



60s 

X+X- 



16s 2 

xix- 
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(B.4) 



T + B = (! + (^ c o + ci/ sg + c 3 / gm + c 4 / sm + c 6 l sq l srn + c 8 l 2 sq 



where 
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C Two hard photon large-angles emission by the initial leptons 



Cross section of 2 photon emission by the initial leptons masses 



e+ (P+) + e (P-) -> 7(Pi) + 7(P2) + hadr(g) 
has a form 



(Cl) 



da 2 ^ 1 a 4 gpp^ dVd^g 
dr h 2!2tt 2 s (g 2 ) 2 wi w 2 



(C.2) 



where the factor 1/2! takes into account the identity of final-state hard pho- 
tons. The relevant contribution to lepton tensor is 



Q 



(2) _ tr 



PPl 



Trp+OjpP-O^p] 



P- 



d 



P1-P2 f^P- 



-12 



d. 



-p + + p 2 „ (j -P+ + Pi „\ ~P+ +P1+P2 

+ [l v ~ A 7 +7 A 1') A Ip 

a_i_2 u+i / w+12 



(C.3) 
(C.4) 
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a_ia + 2 

+ , \ Y'(-P+ + Pihpip- - P2h a , 

a_2d+i 

and 

d_i2 = (p_ -pi -P2) 2 -m 2 ; (C.5) 

d-i = (p~ - Pi) 2 - m 2 ; gL 2 = (p- - P2) 2 - m 2 ; 
d+12 = (~P+ +Pi+ P2) 2 - rn 2 ; 

d+i = (~P+ + Pi) 2 ~ m2 ', d+2 = (p+ + P2) 2 - m 2 . 

Tensor Q^l obey the gauge invariance Q p 2 p \q P = Qp^Qp! = and can be put 
on the form 

Qfl = A 99 PP1 + [A-V-P- + + Anhh + + P-P+jjC.6) 

+A +1 (p + p 1 +pip+) + 

coefficients can be obtained in the standard way: constructing the values 

B g , B n , B ++ , B — , ... = Q P p 1 [g P p 1 , Pi P Pi Pl , p +p p +pi , ...] 
and solving the set of 7 linear equations. 
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